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ABSTRACT 

In an attempt to place an explicit constraint on dark matter models, we define and estimate the mean 
surface density of a dark halo within the radius of the maximum circular velocity, which is derivable 
for various galaxies with any dark matter density profiles. We find that this surface density is generally 
constant across a wide range of maximum circular velocities of ~ 10-400 km s -1 , irrespective of the 
different density distributions in each of the galaxies. This common surface density at high halo-mass 
scales is found to be naturally reproduced by both cold and warm dark matter (CDM and WDM, 
respectively) models, even without employing any fitting procedures. However, the common surface 
density at dwarf-galaxy scales, for which we have derived from the Milky Way and Andromeda dwarf 
satellites, is reproduced only in a massive range of WDM particle masses, whereas CDM provides 
a reasonable agreement with the observed constancy. This is due to the striking difference between 
mass-concentration relations for CDM and WDM halos at low halo-mass scales. In order to explain 
the universal surface density of dwarf-galaxy scales in WDM models, we suggest that WDM particles 
need to be heavier than 3 keV. 

Subject headings: dark matter - galaxies: dwarf spheroidal galaxies - galaxies: kinematics and dy¬ 
namics - Local Group 


1. INTRODUCTION 

Cold dark matter (CDM) theory has serious issues 
on galactic and sub-galactic mass scales, even though 
it has been successful in interpreting a multitude of ob¬ 
servational phenomena at cosmological scales. These in¬ 
clude the so-called missing satellite, core/cusp discrep¬ 
ancy, and too-big-to-fail (TBTF) problems. One of the 
possible mechanisms to resolve these problems within the 
CDM paradigm is to rely on the effects of baryonic pro¬ 
cess such as the photoionization of neutral gas with UV 
background and the impact of gas outflow mainly due to 
mainly supernova feedback on dark halo structures. 

More radically, a possible solution that is well moti¬ 
vated from particle physics is to replace CDM with warm 
dark matter (WDM). The main difference between CDM 
and WDM is that while CDM particles have negligible 
thermal velocities, WDM ones have significant thermal 
velocities during the era of structure formation and thus 
these particles have a free-streaming length, giving rise to 
a cutoff matter power spectrum at larger wave numbers 
corresponding to small spatial scales. The possi ble WDM 
cand idates are the sterile neutrino (see, e. g.. iKusenkc 
2009, for a review) or the gravitino (e.g., iMoroi et al. 
19931 1. These candidates would have a mass in a few keV, 
leading to the suppression of structures on dwarf-galaxy 
mass scales. Therefore, since the cutoff wavelength de¬ 
pends on the particle mass, investigating the nature of 
dark halo structures at low-mass scales leads to useful 
limits on the WDM particle mass, towdm- 

One of the powerful tools for constraining the 
WDM particle mass is Lya forest observa tions (e.g., 
iNaravanan et al.ll2000t IViel et al.ll2008l . 120131) . These ob¬ 
servations allow us to study a matter power spectrum 
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down to small s cales and over a large redshift range. 
IViel et all (12013? ) have observed the flux power spec¬ 
trum of Lyct using the high-resolution spectra and have 
compared it to predicted ones from high-resolution N- 
body and hydrodynamical simulations, assuming that 
WDM masses are 1, 2, and 4 keV. They concluded 
that in order to reproduce the observed power spectrum, 
WDM particles are required to have a lower limit of 
mwDM >3.3 keV (2er confidence level). 

Dwarf spheroidal (dSph) galaxies also provide useful 
tests to distinguish between CDM and WDM or set con¬ 
straints on particle mass of WDM. First, based on the 
difference of dark matter concentration between CDM 
and WDM subhalos, reproducing the maximum circular 
velocity of the Milky Way (MW) satellites can distin¬ 
guish dark matter scenar ios. T his tes t is linked to the 
TBTF problem in CDM. iLovell et all (|2012D performed 
A r -body simulations with WDM and presented that the 
TBTF probl em does not exist wit hi n the frameworks of 
WDM. Also, iPoliscnskv fc Ricottil (120141) showed that a 
WDM particle mass less than 4 keV can mitigate the 
TBTF problem in high-resolution simulations. The sec¬ 
ond test is the method to resolve the missing satellite 
problem. Based on the difference of the abundance 
of sub h alos with in a host CDM and WDM dark halo, 
ILovell et al. (|2014f ) found a conservative lower limit of 
itiwdm >1.1 keV in each V-body simulation. However, 
these tests largely rely on the mass of the MW halo, 
where the estimates of the MW halo mass is still very un¬ 
certain with mass ranging from ~ 0.8 x 10 l2 M^ to ~ 2 x 
10 12 M n (e.g.. [Sakamoto et. al.ll2f)f)3t iDeason et al.l 120121: 

Bovlan-Kolchin et al.ll2013n Kaflc et al.ll20l4 iPifH et Till 

20141) . For the case of the TBTF problem, CDM models 
can resolve the issue i f the mass of the MW halo is less 
than ~ 1.5 x 1O 12 M 0 dWang et al.ll2012l) . For the abun¬ 
dance of subhalos, if the MW halo mass is heavier than 
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1.8 x 10 12 Mq, only mas ses of WDM p articles greater 
than 2 keV are accepted (iKenncdv et al.ll20l4l) . Because 
of the dependence on the assumed MW halo mass and 
the uncertainty in it, the two tests above using the dwarf 
satellites are not yet sufficient to confine CDM or WDM 
scenarios. 

In this Letter, as yet another test on this issue from 
dwarf-galaxy scales, we revisit and consider the dark 
matter surface density, which appears nearly constant 
and independent of galaxy luminosity discovered by 
previ o us studies (e.g.. IDonato et al.l 120091 IGentile et al.1 
20091 ISalucci et all 20121 iKormendv fe Freemanl 12014 
Burkertl 120151) . We propose that the re-definition and 
adoption of this surface density estimated within a ra¬ 
dius of maximum circular velocity directly allow us to 
obtain the limits on particle masses of WDM. This 
work is in part based on the kinematical analysis 
of the MW and Andromeda (M31) dSphs, for which 
more detailed description s will be presented elsewhere 
(IK. Havashi fe M. Chiba 1120151 in preparation). 


2. DARK MATTER SURFACE DENSITY WITHIN A RADIUS 
OF MAXIMUM CIRCULAR VELOCITY 

Recently, analyzing the data of velocity dispersions, 
rotation curves of HI gas, and weak lensing for samples 
of dSph, dwarf irre gula r, spiral, and elliptic al galaxies, 
IDonato et al.l (|2009H and lGentile et al.1 (1200911 fitted these 
data based on cored dark matter density profiles, such as 
the Burkert and pseudo-isothermal profiles, and found 
that the product of the central core density and core 
radius, poAn which is proportional to the central dark 
matter surface density, is constant for all these galaxies. 

While these works offer us important suggestions on 
the universality of the dark halo structure, we note that 
this definition of po r o cannot be applied directly to a 
system with a cusped density profile. Both po and tq 
in such a case may be defined in a specific way depend¬ 
ing on an assumed density profile, so there exists some 
ambiguities when we compare surface densities between 
different systems and also compare with predictions from 
CDM or WDM. 

To adequately compare dark matter distributions be¬ 
ing independent of their density profiles, we introduce a 
mean surface density of a dark halo within a radius of 
maximum circular velocity, 


= 


M(r max ) 


(1) 


where 


M(r max ) = 


4np dul (r')r' 2 dr', 


Pdm(' r ) indicates any dark matter density profiles, and 
r ma x is a radius at maximum circular velocity, K ma x, of 
assumed dark halo profiles, pdm- The surface density de¬ 
fined in equation ([!]) is exactly proportional to the central 
dark matter surface density, Ey max oc p*r*, where p* and 
r* are the central density and scale length of arbitrary 
dark halo profiles, respectively. We note that r max is a 
better-defined scale than a core radius, ro, because the 
latter is uncertain for a cuspy density profile: ro may be 
defined as a half-mass radius or some fraction of it in 
such a case, but it then suffers from the ambiguity of at 



Fig. 1.— Mean surface density of a dark halo, Sy max , within a 
radius at maximum circular velocity, Vmax, as a function of Vkiax 
for different galaxies and Hubble types. While the inset shows all of 
the samples as well as CDM ( WMAP7) prediction, the main panel 
highlights the plots only at Vmax 100 km s"" 1 . Diamonds denote 
the MW and M31 dSphs derived in IK. Havashi &; M. Chiba I (120151 . 
in preparation), where magenta and cyan ones correspond to dSphs 
having a cored and cusped da rk halo profile , respe ctively. Blue and 
red triangles are the original ISpano et al . (2008) sample of spira l 
galaxy data and nearby spirals in THINGS |dc Blok et al. 2008), 
respectively. The data for dwarf irregular galaxies are indicated 
by green circles, and spirals and ellipticals investigated by weak 
lensing are labeled by black squares. Solid lines show the results of 
CDM models with different cosmological parameters of WMAP 7 
(dark green) and Planck (orange) measurements. Red and green 
dashed lines denote WDM models with particle masses of 2 and 
3 keV, respectively. 

what outer boundary the halo is tidally limited, e.g., the 
Navarro-Frenk- White (NFW) profile having a diverging 
total mass. 

Using this definition ©> we evaluate Sy max within 
max taken from the literature based on the HI gas ro¬ 
tation curve of late- an d early-type spirals with pseudo- 
isoth ermal dark halos (jdc Blok et al.(l2008l ISpano et al.l 
I200g), and dwar f irre g ulars with Burkert dark ha¬ 
los ( Gentile et al.1120051 120071) . and the galaxy-galaxy 
weak lensing from spiral and elli ptical g a laxies with 
Burk ert profiles (analyzed by IDonato et al.l 120091 : data 
from iHoekstra et al.l 120051) . For smaller, dwarf-galaxy 
scales, we analyze the stellar kinematics in the seven 
MW (Carina, Fornax, Sculptor, Sextans, Draco, Leo I, 
and Leo II) and five M31 (And I, And II, And III, 
And V, and And VII) dSphs and construct axisymmet- 
ric mass models for these dark-matter-d ominated dwarf 
galaxies (IK. Havashi fc M. Chiba 1 120151 in preparation, 
which con siders more generalized models than our previ¬ 
ous work. (Havashi &; Chiball2012l ). For all samples, ex¬ 
cept for dSphs, we use the dark halo parameters (central 
density and core radius) estimated by the above original 
papers. 

For dSp hs, we employ the axisy mmetric Jeans equa¬ 
tions as in lHavashi fc Chibal (I2012D . but we consider the 
effect of an additional velocity anisotropy of tracer stars, 
(3 Z = 1 — to determine the global shape of a dark 

halo, where vl and v 2 R are vel ocity dispersions of stars 
i n cy lindrical coordinates (See IK. Havashi fc M. Chiba I 
120151 in preparation, for more details). Together with 
[3 Z , we adopt an inner slope of the dark halo density pro¬ 
file as a free parameter as well as the parameterization of 
central density and scale length of a dark halo, apply our 
models to the two-dimensional line-of-sight velocity dis- 
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persion map for the above dSphs, and obtain the best-fit 
dark halo parameters. We have found that whilst most 
of the dSphs have cored or shallower cusped dark halos, 
the Draco, Leo I, And III, and And V dSphs show a steep 
inner density slope; not all of the dSphs have cored den¬ 
sity profiles in th e ir dar k halos, as already suggested in 
lHavashi fe Chibal (1201211 . 

The symbols with error bars in Figure [T] show the es¬ 
timated £y max of the above data sample as a function of 
V max . It is found that even though dark halos in each 
galaxy sample are assumed to have different mass pro¬ 
files (cusp or core) and are estimated by independent 
methods, the dark matter surface density within r max is 
sufficiently constant across a wide range of galaxy masses. 

We discuss probable systematic uncertainties for the 
disk and dSph galaxies, respectively. For the disk galax¬ 
ies, we investigate the impact of the uncertainties in the 
adopted stellar mas s-to-light ratio, M*/L*, on Ey max . 
Ide Blok et al.1 (120081) assumed the three M*/L*s and fit¬ 
ted dark halo parameters for each assumption. Adopt¬ 
ing these parameters, we confirm that these uncertain¬ 
ties in M*/L* do not largely affect the estimation of 
Ey max , varying it by a factor of three at maximum. For 
dSphs, we consider the influences of baryonic components 
and the velocity anisotropy of stars on the estimation of 
Sy max . Since the fraction of the stellar mass in the to¬ 
tal mass of the dSphs is so small, ranging from 10“ 1 
to 10 —3 , that the estimated values of Ey max may not at 
least change beyond a factor of 10, even if these small 
effects of baryonic components are taken into account. 
In our models, we assume a constant velocity anisotropy 
of their member stars and find that this anisotropy is 
almost independent of other model parameters, such as 
a central density and scale len gt h as well as an inner 
slope of a dark halo (IK. Havashi fe M. Chiba 1 12015L in 
preparation). Therefore, the observational evidence for 
the constancy of Ey max is robust. 

3. COMPARISON WITH DARK MATTER MODELS 


The use of the dark matter surface density defined 
by using U max and r max allows us to compare the re¬ 
sults from the data sample with those from cosmological 
A-body simulations within the A-dominated CDM and 
WDM models. We calculate Ey max from the outcomes of 
the A-body simulations as follows. 

First, we adopt the NFW dark matter density profile: 


p{r) 


_ Ps _ 

(r/r s )(l + r/r s ) 2 ’ 


( 2 ) 


which can reproduce CDM (e.g.. iNavarro et al.1 11996L 
Il997f) and WDM fe.g.. lSchneiderl 2012tlLovcll ct al .120141 ) 
dark halos within a wide halo mass range. 

Second, we calculate the scale density, p s , and virial 
radius, r v i r , for a NFW dark halo at the present time 
(z = 0) from cosmological context: 


_ Pcrit A m C 3 

3 ln(l + c) — c/(l + c) ’ 
( 3 M vir ^1/3 

rvir V47rp crit A ro ; ’ 


( 3 ) 

( 4 ) 


where A m = 97 i s a s o-called virial overdensity at 
z = 0 dPrada et al.l I2012h . c is a concentration param¬ 
eter defined as c = r v i r / r s . p cr ; t = 3Hq/8ttG is a critical 



Fig. 2.— Mass-concentration relations for dark halos of CDM 
(orange) and WDM with particle masses of 2 keV (red) and 3 keV 
(green). 


density of the current universe, and Hq denotes the Hub¬ 
ble constant at z = 0. 

Third, using the mass-concentration relation (MCR) 
for a dark halo, we obtain the scale radius, r s . For CDM 
models, we adopt two different MCRs at z = 0 based on 
A-body simulations with slightly different cosmological 
parameters: 


cw = 10.5 ( 


M v 


10 


-0.08 r 


1 + 0 


,15( 


M v 


0.5-i 


c P = 10 


• 7 ( 


M v 


x-O.iir / 

) H 


10 12 h~ 1 Mg 


lO^h^Mo 
(for WMAP7) 
M vir 


, 0.4 


,2.4 x 1O 16 /i _1 M 0 / 


(for Planck) 


where h is a dimensionless Hubble constant. The 
MCR with cosmolog ical parameters of WMAP7 is de- 
riv ed by lPrada c t al.l (12012D . whilst another one is found 
bv IKlvpin et al.l ~ 20141 ). For WDM models, we adopt 
the metho d to derive MCR using the power spectra de¬ 
veloped bv iSchneided (12014 hereafter SI4). In order to 
evaluate the power spectrum of WDM models, the trans- 
fer function formula of the WDM universe bv lYiel et al.l 
(120051) is adopted: 


Twdm(^) = [1 + (ak) 2 ^] 5/,/i , 
with p = 1.12 and 

a = 0.049 


'mwDM' 

-1.11 

Mwdm' 

0.11 

' h ' 

. keV . 


0.25 


L0.7J 


1.22 


( 5 ) 

h -1 Mpc, 

( 6 ) 

where k is a wavenumber and O^dm denotes a cosmolog¬ 
ical density parameter for WDM. Then, the linear power 
spectrum in the case of WDM is derived by 


AvDm(^) = ?W DM (fc)PcDM(fc)- 


( 7 ) 


Using the above power spectrum an d MCR’s fittin g func¬ 
tion described by equation (24) in IKlvpin ct al.1 (12014 , 
MCRs in the case of 2 and 3 keV particle masss of WDM 
are derived by comparing with the A-body simulations 
(see S14 for more details). Figure [2] displays the MCRs 
derived from CDM models with Planck data and WDM 
models with 2 and 3 keV particle masses at redshift 0. 
While CDM models predict that a concentration param- 




























































4 


HAYASHI & CHIBA. 


dSphs (core) [This work] 
dSphs (cusp) [This work] 

CDM (Planck) 
WDM 2keV • 
3keV 



30 

^max [km/s] 


Fig. 3.— Same as Figurebut for dSphs only. Black and gray 
shaded regions are 3cr and 5cr error areas, respectively, for each 
sample. 


eter monotonically increases with decreasing halo mass, 
WDM halos have a clear turnover at a low-mass halo. 

Finally, using the dark matter density profiles obtained 
by the above procedure, we compute I4iax and r max , and 
then calculate Ey max for CDM (WMAP and Planck) and 
WDM (particle masses of 2 and 3 keV) models. The 
solid and dashed lines in Figure |T] are predicted Ey- max 
versus I4iax for four cases of the dark matter models. 
We emphasize that, surprisingly, these theoretical lines 
can well reproduce the data at high-mass scales, even 
though we do not perform any fitting to the data; actu¬ 
ally, the data points (including low-mass scales) appear 
slightly above even the theoretical line giving the high¬ 
est Ey max (CDM with Planck parameters), which may 
suggest there are still some uncertainties in structural 
parameters for dark halos or even cosmology, but the fol¬ 
lowing results for the limit of WDM remain unaffected. 
For all mass ranges, CDM scenarios give a reasonable 
match to the data, where the results based on WMAP 
and Planck data are not so different at low-mass scales. 
This is only due to the slight effect of the difference of 
cosmological parameters on the concentration of dark ha¬ 
los. In particular, V ma . x or r max can depend on both the 
rms amplitude of linear mass fluctuation in a sphere of 
8/i —1 Mpc, a s, and the spectra l inde x of t he initial power 
spectrum, n s dPolisenskv fc Ricottill2014l l. On the other 
hand, WDM ones are inconsistent with the universality 
of Ey max at dwarf-galaxy mass scales. This different fea¬ 
ture between CDM and WDM models at the low-mass 
end stems from the properties of the halo concentrations 
at low halo-mass scales. From Figure [2] the concentra¬ 
tions of WDM halos differ substantially from those of 
CDM halos at low halo-mass scales, thereby the value of 
Ey max decreases prominently at small I4iax correspond¬ 
ing to dwarf galaxies. 

We also investigate how far WDM models deviate from 
the data. Figure [3] is the same as Figure [Q but zooms in 
the dwarf-galaxy scale and plots 3 a and 5cr error areas of 
each sample. Focusing on I4iax < 20 km s -1 where the 


difference in Ey max emerges significantly, even a 3 keV 
particle mass of WDM appears to deviate beyond a 30- 
error in our sample. Following this result, in order to 
be in agreement with the universality of Ey max from the 
data, WDM particles need to be heavier than 3 keV. 

4. SUMMARY AND CONCLUSIONS 

In this Letter, we define and estimate a mean surface 
density of a dark halo within a radius of maximum circu¬ 
lar velocity, Ey max and apply it to dark halos in various 
galaxies. This physical value reflects the central dark 
matter surface density, thus not strongly affected by the 
influence of external disturbances. We find that even 
though dark halo parameters in each galaxy sample are 
derived based on the assumption of different dark halo 
density profiles and determined by several independent 
methods, this surface density is universal across a wide 
luminosity range of galaxies, irrespective of dark matter 
distribution in each of the galaxies. 

Furthermore, this constancy is found to be useful for 
constraining dark matter models, especially the parti¬ 
cle mass of WDM, because Ey max is directly derivable 
from the theoretical predictions based on iV-body simu¬ 
lations. Our results show that at larger values of y rnax 
corresponding to high halo-mass scales, this constancy 
for real galaxies can be successfully reproduced in both 
CDM and WDM models, even though we do not per¬ 
form any fitting to the data. However, Ey max derived 
from WDM largely deviates from the observational con¬ 
stancy at dwarf-galaxy mass scales, whilst that of CDM 
is in good agreement with it. This is because there is a 
prominent difference between MCR for CDM and WDM 
at low halo-mass scales. In order to be in agreement with 
this constancy, we have found that WDM particles need 
to be heavier than 3 keV. 

This Letter proposes an astrophysical test on the na¬ 
ture of dark matter. The key dark halo mass range to 
distinguish dark matter models in this test is the lowest 
halo-mass scales, such as those of classical and ultra-faint 
dwarf galaxies. Wide-held spectroscopic studies of these 
galaxies over large areas with the planned Prime Focus 
Spectrograph to be mounted on the Subaru Telescope 
and high-precision spectroscopy for a plenty of much 
fainter stars with the Thirty Meter Telescope will pro¬ 
vide a better determination of dark halo properties of 
the MW and M31 dSphs, and thus allow us to increase 
our knowledge about the basic properties of dark matter. 
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